Effect of rolling on dissipation in fault gouge 
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Sliding and rolling are two outstanding deformation modes in granular media. The first one 
induces frictional dissipation whereas the latter one involves deformation with negligible resistance. 
Using numerical simulations on two-dimensional shear cells, we investigate the effect of the grain 
rotation on the energy dissipation and the strength of granular materials under quasistatic shear 
deformation. Rolling and sliding are quantified in terms of the so-called Cosserat rotations. The 
observed spontaneous formation of vorticity cells and clusters of rotating bearings may provide an 
explanation for the long standing heat flow paradox of earthquake dynamics. 
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I. INTRODUCTION 

The micromechanics of heat production by friction in 
granular materials has become an important issue in the 
study of earthquake mechanics. One of the unresolved 
controversies in this field is a phenomenon that geophysi- 
cists call the heat-flovi' paradox Jj . According to common 
sense, when two blocks grind against one another, there 
should be friction, and that should produce heat. How- 
ever, measurements of heat How during earthquakes are 
unable to detect the amount of heat predicted by simple 
frictional models. Calculations using the value of rock 
friction measured in the laboratory, i.e. a typical friction 
coefficient between 0.6 and 0.9 0, lead to overestimation 
of the heat flux. As an example one refers in this con- 
text to the heat flow observations made around the San 
Andreas fault, which show that the effective friction co- 
efficient must be around 0.2 or even less One possible 
scenario for the explanation of these observations is the 
mechanism of heat-induced pore-fluid pressure increase 
. Other mechanisms have been also discussed 
At any rate, with or without pressurization, the 
correct assessment of frictional heat production during 
shear is a central issue. Here we will address this issue 
by resorting to the micromechanics of dry granular media 
representing the gouge, i. e. the shear band consisting of 
fragmented rock inside the fault zone. 

The formation of rolling bearings inside the gouge has 
been introduced as a possible explanation for a substan- 
tial reduction of the effective friction coefficient [a . This 
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simplified picture assumes that the gouge is filled with 
more of less round grains which, as the plates move, can 
roll on each other thus reducing the amount of frictional 
dissipation. Granular dynamics simulations Q, |3j E)] ^nd 
Couette experiments [l^ have demonstrated the sponta- 
neous formation of such bearings. 

The overall effects of grain rotation are studied here 
using granular dynamics simulations. We show that par- 
ticle rotation induces a phase separation in the granular 
media in terms of three coexisting phases: (1) Vorticity 
cells, (2) rotational bearings and (3) slip bands. The first 
two phases reduce significantly the frictional strength and 
the dissipation with respect to the hypothetical case of 
simple shear. We quantify these phases in terms of the 
so-called Cosserat rotations. We address the necessity to 
introduce these variables for the constitutive modeling of 
fault gouge. In Sees Ull and IHII we present the theoreti- 
cal background and the particle-based model. In Sec. IIVI 
the effect of particle rotation on the strength and fric- 
tional dissipation is investigated. In Sec. we calculate 
the population of the three coexisting phases using the 
homothetic-antithctic decomposition of the Cosserat ro- 
tations. 



II. COSSERAT CONTINUUM 

In the framework of continuum mechanics, the math- 
ematical description of granular rolling and sliding is a 
challenging task. We notice first that classical contin- 
uum theories introduce the concept of the material point 
as a representative assembly of grains and ascribes to it 
only the degrees of freedom of displacement, which in 
turn are correlated to the displacements of the grains of 
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the assembly. Thus classical continuum theory makes 
no provision for particle rotation. More recent contin- 
uum models include the rotational degrees of freedo m by 
using the so-called Cosserat rotations 0, 0, 0, IT^ . 
These are continuum field variables, measuring the rela- 
tive particle rotations with respect to the rotational part 
of the displacement-gradient field. The name of these 
variables gives tribute to the brothers Cosserat (1909) 
who were the first to propose such a continuum the- 
ory. 50 Years after the first publication of the original 
work [10], the basic kinematics and static concepts of 
Cosserat continuum were reworked in a milestone paper 
by Giinther Giinther's paper marks the rebirth of 

micromechanics in the 1960s. Following this publication, 
several hundred of papers were published on the subject 
of micromechanics of granular media _16\ . The growing 
interest in the Cosserat theories in recent years followed 
the link that was made by Miihlhaus and Vardoulakis 
[T^ between the Cosserat continuum and the onset of 
shear-bands in granular materials. Observations of par- 
ticle rotations from particle-based and continuum-based 
computer models |l7l | and physical experiments jlSj show 
that inside shear- and interface-bands the particles rotate 
differently as their neighborhood. These findings demon- 
strated the necessity to introduce the Cosserat rotations 
as additional field variables in the shear-band evolution. 
These new variables involved also a characteristic ma- 
terial length that allows in turn to reproduce the char- 
acteristic width of shear bands. This internal length has 
special significance from the computational point of view, 
because its resolves the mesh- dependency problems in the 
Finite Element simulations [l^l20j|. 

At any material point of the Cosserat continuum we 
assign a velocity v and a spin vector uj . Accordingly 
in plane-strain deformation the kinematic fields are: the 
classical strain-rate tensor, which corresponds to the 
symmetric part of the average particle velocity : 

and the Cosserat rotation, that is given by the difference 
between the macro- and micro-rotations: 



tions. The discrete model consists of disks confined be- 
tween two horizontal plates. A normal force is applied on 
the plates, as they are moved in opposite directions with 
a constant velocity. Periodic boundary conditions are 
imposed along the horizontal direction, see Fig. ^ This 
geometry is a simplified model of a gouge. Real rock 
gouge consists of non-spherical particles. In the present 
model we are able to simulate two extreme cases: The 
first one represents a young fault, which is characterized 
by a strong interlocking of closely packed rocks in the 
gouge. This case is simulated by hindering the rotation 
of the disks. In the second case we allow the particles to 
rotate without rolling resistance. This is an idealization 
of mature gouges, where the interlocking and the rock 
asperities are reduced due to grain fragmentation. 

The discrete model is a 2D implementation of the Lat- 
tice Solid Model Hllllil. This is a suitable platform 
to investigate the dynamics of fault gouge via granular 
dynamics simulations. Two disks of radii Ri and Rj in- 
teract when the distance between their centers of mass 
rij is less that the sum of their radii. Their interaction 
is given by the viscoelastic contact force 



= knAxnfi + ktAxtt + Som{-/nVnn + ItVft) (3) 

The first two terms are elastic forces, and the last two 
are viscous forces. The unit normal vector n points in 
the direction of the vector connecting the center of mass 
of the two disks. The tangential vector t is taken per- 
pendicular to n. The elastic material constants are the 
normal fc„ and tangential kt grain contact stiffnesses, 
The normal elastic deformation is the overlapping length 
Aa;„ = Ri + Rj — rij . The tangential elastic deformation 
is chosen to be consistent with the Coulomb sliding con- 
dition as follows: When two particles come into contact 
we set Axt — 0. Then, at each time t, we guess a new 
value for the tangential elastic deformation as: 



AxP {t) = Axt {t - dt) + v^tdt, (4) 



l,dux duy 
^=2(^-9^)-"' 

This variable can be calculated as the difference between 
the rotation of the branch vector - defined as the line 
connecting particle mass centers - and the rotation of 
the particle, averaged over all pair contacts in the repre- 
sentative volume element. 



III. PARTICLE-BASED MODEL 

We will investigate the discrete counterpart of the 
Cosserat rotations by using granular dynamics simula- 




fPo 

FIG. 1: Contact network in the shear cell at critical state. 
The lines represent the branch vectors; the width of the lines 
the normal contact force. 
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where is the tangential relative velocity at the contact: 

= LUiRi + ojjRj + [vi - Vj) ■ t. (5) 

Here Vi is the velocity and LOi is the angular velocity of 
the particles in contact. The predicted value of elastic 
deformation should be corrected to satisfy the Coulomb 
sliding condition |/f | < where and are the 

tangential and normal component of the elastic force and 
fj, is the friction coefhcient. This condition equates to: 



of the plates vq = 10~^. The time step is At = 0.05. 
The most relevant parameter of this model is the ratio 
between the shear velocity vq and the velocity of com- 
pressional waves, which in or model is Vp ^ Ro\/kn/m. 
In our simulation ~ 1 so that Vp/vQ ^ 10^. This value 
should be compared to the ratio in realistic fault zones, 
where Vp ~ IKm/s and vq ^ Icm/year leading to a fac- 
tor of Vp/vo ~ 10^". To remedy this discrepancy of time 
scales we use the quasistatic limit: The velocity is chosen 
low enough so that the reduction of it by half affects the 
effective friction coefhcient by less than 5%. 



Axt{t) = sign{AxP{t))mm C " , ' JAxfl). (6) 

kt 

The viscous force in Eq. Q assures the restitution be- 
tween colliding particles. 7„ and jt are coefficients of 
viscosity, and the harmonic mean m — {l/mi + l/mj)^^ 
is the effective mass of the disks. The normal and tan- 
gential components of the relative velocity at the con- 
tacts are = {vi — Vj) • n and given by Eq. (jSJ. The 
factor 5o = Ax„/Ro is included in Eq. 13 to guarantee 
continuity of the viscous force during collision. Rq is the 
averaged radius of the disks. 

The roughness of the driving plates is modeled by at- 
taching particles to it. Their vertical displacement is 
controlled with a simple viscoelastic force — knS + 
"fnTriiV^^ where S is the overlapping length and is the 
relative normal velocity. The attached disks are not al- 
lowed to rotate and their horizontal velocity is set to the 
velocity of the plates. 

Each contact contributes to a direct force and a 
torque r"^ — Rif^-t in the equation of motions. The 
model does not include gravitational forces, but a viscous 
force /" = 'yrriiVi and a torque r = ^miRfuji is included 
for each particle. This viscous forces allow relaxation of 
particles without contacts. We use the Verlet method for 
solving the equations of motion ^2^. 

The efficiency of the simulation is mainly determined 
by the method of contact detection. Our method searches 
in each step the contacts in a list of neighbours that is 
called a Verlet List. This list is constructed by taking 
the pair particles whose distance of their center of mass 
satisfies the constraint r'^ < Ri + Rj + S. The Verlet 
List is updated when the maximal displacement of the 
particles since the last update is larger than 6/2. The 
parameter 6 is chosen by making a compromise between 
the storage (size of the Vertex List) and the compute time 
(frequency of list updates). A Linked Cell algorithm is 
used to allow a rapid calculation of the new Verlet List 

M 

The material constants of the model are the normal 
stiffness kn — 1; the tangential stiffness kt = 0.1; the 
normal 7„ = 0.001 and tangential 7„ — 0.0001 damp- 
ing frequency and the body viscosity 7 = 0.0001. The 
density of the disks p = 1; and the friction coefhcient 
whose default value is fi = 0.5. The control parameters 
are the applied pressure po ~ 0.001, and the velocity 



IV. EFFECT OF PARTICLE ROTATION 

Here we address the question of how particle rotation 
affects the mechanical response of the shear cell. Simu- 
lation of shear cells with rotating and non-rotating disks 
arc compared by calculating the power dissipation and 
the stress ratio at the limit of large shear deformation. 
Each shear cell consists of 1600 disks with random sizes 
between 0.5i?o and 1.5i?o. The length of the cells is 80i?o. 
We start from a loose packing that is compressed by ap- 
plying a constant pressure at the top plate. After a short 
collisional regime the sample reaches an equilibrium con- 
figuration where the kinetic energy decreases exponen- 
tially with time. After this stage the sample is sheared by 
applying a horizontal velocity vq at the top plate and —vq 
at the bottom. The simulations are performed by tak- 
ing microscopic friction coefficients between 0.0001 and 
8. These values should be compared with the friction of 
realistic materials, which ranges from 0.0001 for super- 
lubricated surfaces to 1.2 for rubber-concrete contact sur- 
faces. 

All samples reach a limit state for large shear deforma- 
tion. This state resembles the so-called critical state of 
soil mechanics As shown the Fig.|21 the stress and 
the void ratio reach a constant value besides some fluctu- 
ations. These fluctuations result from the characteristic 
stick-slip dynamics of the shear cell: In the stick phase, 
the elastic energy is stored in force chains. These chains 
build up during the slow relative motion of the plates. 
The elastic energy is released in form of quakes. Each 
quake corresponds to the collapse of a force chain which 
is reflected by an abrupt drop of the macroscopic stress 
and a sudden compaction of the sample. The collapse of 
force chains leads to reorganization of the particles and 
hence acoustic emission. This is detected from the accel- 
eration of a single particle in the sample, see part (d) of 
Fig. 121 Between two quakes the buckling of forces chains 
leads to an overall elastoplastic response with constant 
increase of void ratio. 

The strength of the shear cell in the critical state can be 
quantified by an effective friction coefficient /i* = Ft/ Fn. 
Here Fn and Ft are the normal and tangential force acting 
on the top plate. Part (a) of Fig. |21 shows the time evolu- 
tion of /I*. Typically this value increases in the interval 
between two quakes, and it drops during the quakes, lead- 
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FIG. 2: Time profiles of (a) ratio between the normal and 
tangential force acting on the plates, (b) void ratio, calculated 
as e — {A — Ad)/ Ad, where A is the area of the cell and Ad 
is the total area occupied by disks, (c) frictional dissipation 
normalized by /lo = LpoVotc and (d) acceleration of one disk 
in the center of the cell, normalized by ao — L/t1. The time 
is normalized by = L/vq 



ing to time fluctuations which are of the same order as 
the time average value /is = /i* . 

The dependence of /is on the microscopic friction co- 
efficient for rotating and non-rotating disks is plotted in 
part (a) of Fig. |3| For small values of /i the strength 
of the shear cells is not affected by grain rotation. In 
both cases the effective friction coefficient is larger than 
/i. This is due to the addition of an interlocking strength 
component to the shearing friction component. In the 
limit case /i ^ the effective friction tends to 0.1, show- 
ing that interparticle friction is not the only origin of the 
macroscopic frictional behavior of granular materials. 

For large values of /i particle rotation has a signifi- 
cant effects on the strength of the shear cell. Particu- 
larly, in the range 0.28 < fi < 1.42 samples with rotating 
disks have a macroscopic friction coefficient lower than 
the contact friction coefficient. On the other hand, the 
macroscopic friction coefficient is larger than the contact 
friction coefficient for non-rotating disks. In the limit 
fi ^ oo, the strength of both samples tends to a con- 
stant value, This value is 0.34 for rotating disks samples 
and 1.43 for non-rotating disks. 

We now turn to the effect of grain rotation on dissi- 
pation. Energy loss in the sample results from frictional 
and viscous dissipation. In our simulations, viscous dis- 



sipation takes place only in the coUisional regime, which 
applies only during the phase of sample preparation. In 
the quasistatic phase and as the sample reaches the criti- 
cal state, the frictional dissipation is at least three orders 
of magnitude larger than the viscous dissipation. There- 
fore the produced heat can be calculated as: 

h = J P{i)dt, (7) 

where P is frictional power dissipation. This quantity is 
calculated as the sum of the tangential force times the 
sliding velocity of all the contacts. 




FIG. 3: Dependence of the (a) effective friction coefficient 
and (b) frictional power coefficient on the microscopic friction 
coefficient for rotating and non-rotating particles. Data of 
five different samples are superposed. Solid lines represent 
the linear interpolation around the averaged value, dashed 
line corresponds to /is = /i and /i/i = fi. 
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FIG. 4: Simple shear model of gouge deformation. 



where is the tangential contact force; the tangential 
relative velocity at the contact and Ax^ the elastic part 
of the tangential displacement as defined by Eq. 
The heat versus time is plotted in part (c) of Fig. |21 for 
rotating particles with fi — 0.5. The general trend of 
the heat is to increase slowly between two quakes, and 
rapidly during quakes. 

The calculated heat should be compared to the theo- 
retical energy dissipation of simple shear, which reads: 

hth = fJ-FnvAt, (9) 

where /i is the microscopic friction coefficient, Fn the 
normal force on the plates, v = 2vo the relative velocity 
of the plates and At is the time interval. The value of 
Eq. (PI is calculated according to the model shown in 
Fig. ^ This model corresponds to the simple shear of a 
multilaminate with interlaminate friction equal to fi. To 
compare this value to the calculated power dissipation 
P{t), we introduce the frictional power coefficient as fih — 
P{t)/Po, where Pq = F„u and P{t) is the time average 
along the critical state regime. Note that fif^ = fi when 
the measured heat coincides with the theoretical value. 

Part (b) of FigOshows the dependency of the frictional 
power coefficient on fj,. Heat production assumes a max- 
imal value at /i = 0.5 for rotating particles and /.i — 1.0 
for non-rotating particles. In both cases the dissipation 
tends to decrease for fj, —>■ and — > cx). In the first case 
because ft^O and in the second case because the frac- 
tion of sliding contacts is very small for large microscopic 
friction coefficients. 

For small values of /i heat production is not affected 
by particle rotation, but for larger values it is affected by 
almost one order of magnitude. In both cases heat pro- 
duction is quantitatively different from the value of sim- 
ple shear. Small friction coefficients lead to a dissipation 
larger than the theoretical one. This is due to the large 
population of sliding contacts in this regime, and their 
relative orientation with respect to the principal direc- 
tion of load, which leads to tangential forces larger than 
the expected value for simple shear. For large friction co- 
efficients few contacts can reach the sliding condition, so 



that the dissipation is much lower than that one expected 
for simple shear. 

The displacement field confirms that the picture of sim- 
ple shear (Fig. ^ is also not supported. In part (a) of 
Fig. ]E\ we plot the displacement field for non-rotating 
particles when /i = 0.5. We observe the formation of 
blocks of particles moving approximately as a whole. The 
boundaries between these blocks are not fiat, but curved. 
This produces strong interlocking, that explains the con- 
siderable increase of the strength of the material with 
respect to simple shear. Note that these blocks can not 
slide against each other for very long without changing 
shape. The reorganization of these structures appears by 
means of large quakes. These quakes change completely 
the displacement fields, even when the position of the 
particles stays approximately the same. 

If the particles are allowed to rotate, vorticities ap- 
pear spontaneously. They are shown in part (b) of FigO 
These vorticities are accompanied by strong temporal 
fluctuations of the displacement field, but they can ap- 
pear and disappear intermittently in the same zone. The 
vorticity field has been observed in many numerical simu- 
lations [23,123. It resembles to some extent the turbulent 
regime observed in fluids, but their dynamics is quite dif- 
ferent p^ : Fluids under slow motion present a laminar 
regime where the mass displacements can be considered 
as simple shear. On the other hand, our shear cells de- 
velop large vorticities even in the quasi-static regime, rul- 
ing out the simple shear deformation regime. The short 
lifetime of these vorticities, typically the same as the in- 
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FIG. 5: Snapshot of the displacement field for (a) non- 
rotating and (b) rotating grains. The friction coefficient is 
Ai = 0.5. 
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The latter approximation is valid when Ax„ << d, where 
Axn is the overlapping length and d the characteristic di- 
ameter of the disks. The unit vector z**^ is perpendicular 
to the plane of the disks, and the unit tangential vector 
satisfies = z^^ x n'^ . 

Let us consider two points attached to each particle, 
in a region infinitesimally near to the contact point. The 
tangential velocity of these points are given in terms of 
the linear Vk and angular oj^ velocities of the particles: 



FIG. 6: homothetic (a) and antithetic (b) couple of the 
Cosserat rotations. 



terval between two large quakes, contrasts to the large 
lifetime of the eddies in turbulent flow. Thermal mea- 
surements in dynamic shear-banding on metals also re- 
veals that homogeneous shear is not possible and vortices 
are the rule |29j| . 

The spatial distribution of vorticity cells, when com- 
bined with the distribution of rolling, provides a novel 
picture of bearings: Inside vorticity cells all particles ro- 
tate as a rigid body, whereas the space between the vor- 
ticity cells is characterized by intense relative rotations. 
This leads to clusters of rotational bearings and zones 
of intense slippage. We will quantify the contribution of 
these modes to global deformation by performing a kine- 
matic decomposition of the contact deformation. This 
decomposition will distinguish rolling from sliding and 
from the rigid body motion of the vorticity cells. 



V. MEASURE OF ROLLING 

Several definitions of rolling can be found in the liter- 
ature [33, 1^ ■ They are taken as measures of contact 
deformation for each pair of interacting particles. Kuhn 
and Bagi decompose the degrees or freedom of the two 
particles into rigid body motion and objective motion 
|31| . Then the rolling is defined as the average of the ob- 
jective translation at each side of the contact point [s^l . 
Based on these studies, we introduce a new definition of 
rolling in terms of the rotation and translation of the par- 
ticles in contact. This definition will be consistent with 
the fact that rolling allows deformation in the assembly 
without frictional dissipation or accumulation of elastic 
energy. 

Let us suppose that at time t two particles indexed by 
i and j are in contact, and they stay in contact during an 
infinitesimal time interval dt afterward. We introduce a 
system of coordinates n^^ , t^^ , z^^ attached to the contact. 
The unit normal vector n^^ connects the center of mass 
of the particle i to that of the particle j: 



Xj (t) - Xi {t) Xj{t) - Xi (t) 



\Xj{t) - Xi{t)\ 



Ri ~\- Rj 



(10) 



Vii- = Vi-t + ujiRi 
'"It = Vj -t-ujjRj. 



(11) 



Let us take the movement of the branch vector, i.e the 
vector connecting the center of mass of the two particles. 
The tangential component of the velocity of a point at- 
tached to the branch vector at the contact will be called 
rigid body velocity: 



rb 



{viRj 



VjRi 



Ri + Rj 



(12) 



If the two particles move as a rigid body, the velocity 
coincides with Eqs. Ijllll . Otherwise there is a relative 
velocity between the two points attached to the particles. 
This velocity can be calculated by subtracting the rigid 
body like velocity of Eq. (|12|l from the contact velocities 
of Eqs. lfTT|l : 



TO 

rb 



(13) 



These velocities can be given in terms of linear and an- 
gular velocities by using Eqs. (|11() and l|12() : 



R^ 



-LUiR 



Ri + R- 
R. 



{v,~v,)-?^. (14) 



They will be called objective velocities of the two particles 
at the contact. These velocities are objective in the sense 
that their magnitude is unaffected by the common rigid- 
body motion of the particle pair. In particular, these 
measures are not affected by a change of the reference 
frame. The objective velocities vanish if, and only if, the 
two particles move as a rigid body. Otherwise the objec- 
tive velocities involve rolling, sliding or accumulation of 
shear strain at the contact. 

We will relate these objective velocities to the afore- 
mentioned Cosserat-continuum rotations. At the mi- 
cromechanical level, the later are defined as the relative 
rotation of the particle with respect to the rotation of 
the branch vector around the axis parallel to z^^ at the 
contact point js^ : 
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UJi 



(15) 



where 4>i is the angular velocity of the particle i with 
radius Ri. O'^ is the rotational velocity of the branch 
vector: 



Ri + Rj 



(16) 



Replacing Eq. (|16|l into Eqs. (|15|l and comparing the 
result with Eqs. (|14|) we obtain a relation between the 
Cosserat rotations and the objective velocities. 



R,W'- 



'RjW^' 



(17) 



These equations provide an interesting connection be- 
tween the relative orientation of the Cosserat rotation 
and rolling: When the Cosserat rotations are in a ho- 
mothetic couple as shown part (a) of Fig. ^ the objec- 
tive velocities s'-' and s-'' have opposite signs. In this 
case, depending on whether the contact is or not in the 
sliding condition, the contact deformation results either 
in frictional dissipation or in accumulation of elastic en- 
ergy, without rolling deformation, nomothetic couples 
appear in several cases: If the two disks in contact rotate 
in the same sense without linear velocity, the Cosserat 
rotations reads W^^ = LOi and W^^ = uij. If the parti- 
cles do not rotate, the Cosserat rotations are given by 
W^^ — W^"^ — SI'-' . In both cases they are in a homoth- 
etic couple, and hence there is no rolling deformation. 

The appearance of rolling between the particles is 
linked to the antithetic couple of the Cosserat rotations. 
This couple is shown in part (b) of Fig. El In this case 
the rolling is given by the common part of the objec- 



tive velocities, i. e. min(|s'' 



In particular, if the 



Cosserat rotations satisfy WiRi + WjRj — both grains 
have the same objective velocity. Therefore the contact 
deformation does not accumulate elastic deformation or 
produce frictional dissipation, and hence the deformation 
corresponds to pure rolling. 

In the general case the objective velocities should be 
decomposed into rolling and dislocation. The first one 
results in: 



yij i homothetic couple 

roll 1 min(|s*-' I, |s-'*|)sign(s'^) antithetic couple, 

(18) 

and the dislocation is given by the difference of the ob- 
jective velocities: 



(19) 



( , it is easy to prove 



Using the identity s*-' — 

that the dislocation velocity corresponds to the relative 




FIG. 7: Snapshot of (a) vorticity field, (b) objective veloci- 
ties at the contacts as shown the inset in part (a), and their 
decomposition in (c) dislocation and (d) rolling. The sample 
consists of rotating grains with ^ = 0.5. 



tangential velocity defined in Eq. JSJ. In the limit of rigid 
disks, where elastic deformation at the contact is absent, 
the dislocation velocity coincides with the sliding veloc- 
ity. In granular dynamics an elastic shear deformation is 
allowed, so that the sliding velocity is given by: 



slid 



disl 



djAxt) 
dt 



(20) 



where Axt is the elastic part of the tangential displace- 
ment at the contact, that is given by Eq. ©. Accord- 
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ing to Eq. ^ the sliding velocity times the tangential 
force corresponds to the frictional dissipation. The sec- 
ond term of Eq. 1201 involves accumulation of elastic en- 
ergy. 

The Cosserat rotations turn out to be a suitable or- 
der parameter to describe the vorticities and bearings in 
the shear deformation: As shown the Fig. [7| the spa- 
tial distribution of these structures is not just random, 
but strong correlations appear in form of three coexist- 
ing phases: (1) Vorticity cells, where the particles rotate 
almost as a whole, so that the Cosserat rotation is vanish- 
ingly small, (2) clusters of particles with intense rolling 
(rotational bearings), where the Cosserat rotations are 
antithetic, and (3) zones between particles with intense 
dislocation (microbands), characterized by homothctic 
couples of the Cosserat rotations. Microbands induce 
frictional dissipation and accumulation of elastic energy, 
whereas the rotational bearings accommodate the vortic- 
ity cells to make them more compatible with the imposed 
kinematic boundary conditions. 

The microscopic friction coefficient can serve as a con- 
trol parameter of the relative population of these three 
phases. This is shown in Fig. IHlfor (a) non-rotating and 
(b) rotating particles. In both cases the sliding displace- 
ment and rigid body velocities are the dominant defor- 
mation modes for small values of fj,. In this regime, the 
dislocation velocity coincides with the sliding velocity, 
because almost all the contacts are in the sliding condi- 
tion. Deformation in this case is characterized by small 
clusters of particles moving against each other through 
microbands of intense sliding. For rotating particles the 
rolling is much lower than the other deformation modes, 
whereas the rolling is absent for non- rotating grains. 

For large friction coefficients the sliding velocity be- 
comes negligible, because few contacts are able to reach 
the sliding condition. In this regime, depending on 
whether the particles are or not allowed to rotate, the 
population of the deformation modes is quite different: 
for non-rotating particles, the contact deformation is 
dominated by rigid-body motion and elastic dislocation. 
The latter one builds up elastic energy that is liberated in 
the form of strong quakes. The dependence of dislocation 
and rigid-body motion on ^ is fairly weak. This reflects 
a self organization of the shear cells, characterized by a 
non-dependency of the effective friction coefficient on fi, 
already shown in part (a) of Fig. O 

For rotating grains and large friction coefficients rolling 
plays a relevant role, as shown part (b) of Fig. |H1 The 
deformation is dominated by rigid-body motion (due to 
the vorticities), rolling (due to the rotational bearings) 
and elastic dislocation (due to building of force chains). 
The self organization of the shear cells is given by the fact 
that rolling and vorticities are not affected by a change in 
the microscopic friction coefficient. The sliding turns out 
to be much smaller that the other modes, because only 
few contacts can reach the sliding condition. Therefore 
the dislocation results almost completely in accumulation 
of elastic energy that is released in the form of quakes. 
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FIG. 8: Decomposition of kinematic roles for (a) non-rotating 
and (b) rotating particles, as function of fi. They are calcu- 
lated form the average of the absolute value of the rolling, 
dislocation, sliding and rigid body velocities at the contacts. 
Data of five diflterent samples are superposed. 



The details of the dynamics of one quake is shown in 
Fig. |51 We use a friction coefficient oi fi — 0.5. We 
observe a small time interval (shock) characterized by a 
sudden decrease of the effective friction coefficient, along 
with an abrupt compaction and a rapid generation of 
heat. This is followed by a longer time interval (after- 
shock) given by an exponentially decay of microseismic 
activity with almost no frictional dissipation. More de- 
tails of the structure of the quake are visible by calculat- 
ing the time evolution of the contact modes of deforma- 
tion. Part (e) of Fig. shows that most of the contact 
deformation at the shock corresponds to rolling and rigid 
body deformation. In the after-shock the sliding veloc- 
ities become vanishing small, whereas the other modes 
decay exponentially. Acoustic waves are generated at 
the point where the force chain fails. These waves travel 
though the sample and they are reflected as they reach 
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FIG. 9: Time evolution of (a) stress, (b) void ratio, (c) fric- 
tional dissipation, and (d) acoustic emission at the center of 
the sample, (e) Average of the absolute value of rolling (cir- 
cles), dislocation (asterisks), sliding (crosses) and rigid body 
velocities (diamonds). The sample consists of rotating parti- 
cles with fi = 0.5. 



the plates. After many reflections they becomes uncorre- 
lated, which removes almost all contacts from the sliding 
condition. Therefore shding velocity modes are not ac- 
tive in the aftershock, whereas the other modes decay 
exponentially in form of uncorrelated oscillations. 

As far as earthquakes is concerned, it is remarkable 
that only a minute part of the contact deformation dur- 
ing the quake corresponds to sliding. This leads to an 
effective friction coefficient around 0.25 which is lower 
than the contact friction cocfhcicnt oi fi — 0.5, see part 
(b) of Fig. O Therefore the introduction of rolling in the 
gouge dynamics could potentially explain the Heat Flow 
Paradox. 



VI. CONCLUDING REMARKS 

The most important contribution of this work is to 
show that spontaneous formation of vorticity cells as- 
sisted by rotational bearings is the mechanism of reduc- 
tion of strength and frictional dissipation in shear cells. 
For the range of rock friction of ^ = 0.6 — 0.9, the effective 
friction coefficient increases from = 0.3 to /is = 0.35, 
whereas the frictional power coefficient decreases from 
/i/i = 0.28 to /i/i = 0.22. These results could explain the 
heat flow anomalies in fault gouges, where the macro- 
scopic friction coefficient, as deduced from heat flow ob- 
servations, is around one order of magnitude lower than 
the rock friction The existence of such rotational 
patterns confirm earlier speculations about the effect of 
rolling in the reduction of heat production by means of 
a self-organization process in fault gouge 1] . 

The effective friction coefficient and frictional power 



coefficient are strongly dependent on particle rotation. 
This reflects the necessity to introduce rotational degrees 
of freedom in the continuum description of fault gouge. 
As deduced from the deformation field at the contacts, 
the Cosserat continuum approach should be consistent 
with the observed three phase separation of kinematic 
modes: (1) vorticity cells, where the Cosserat rotations 
are absent; (2) Bearings given by antithetic couple of the 
Cosserat rotations, and therefore a pronounced rolling 
at the contacts; and (3) microbands of homothetic cou- 
ples, with pronounced dislocations, and hence, strong ac- 
cumulation of elastic energy and high frictional dissipa- 
tion. This description requires an extension of the ex- 
isting Cosserat continuum models, which consider only 
homothetic couples of the Cosserat rotations 'T?,'2^. As 
far as dissipation is concerned, besides the Cauchy stress 
tensor, a couple stress tensor should be introduced as 
an additional static variable of the enhanced continuum, 
entering as the energetically dual counterpart of the gra- 
dient of the homothetic part of the Cosserat rotations 
[ill IT^ [T3 . Is^ l. In this context numerical simulations 
can be used as a virtual laboratory to assist the develop- 
ment of micromechanical constitutive models. 

An improved constitutive law for fault gouge will re- 
quire integration of not only rolling, but also the con- 
tribution of grain fragmentation to the energy budget. 
Observations of fault zones suggest that earthquakes can 
pulverize rocks in the gouge 01 ■ This leads to a frac- 
tal grain size distribution and an increase of the gouge 
surface area. Some important issues concerning the heat 
flow paradox should be considered: Does fracture energy 
play an important role in the earthquake energy balance? 
and do such fractal gouges develop vortical structures and 
rotating bearings?. 

There are also some aspects about the aseismicity of 
fault zones which deserve detailed study in a future work. 
Vortical structures and rotating bearings promote a co- 
herent deformation which remains during two quakes. 
These may inhibit the large events we observe in samples 
with non-rotating grains. The dynamics of these rota- 
tional patterns will be significant in the understanding 
the enigmatic aseismic creep, where two tectonic plates 
move against each other without accumulating elastic en- 
ergy or generating earthquakes (35. .36.] . 
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